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ABSTRACT 

An e ± atmosphere is inevitably created around a black hole accretion disc, the spectrum of 
which extends to MeV energies. Pairs created in 7 — 7 collisions outside the disc are blown 
away by soft radiation (which dominates the bolometric luminosity of the disc) and form a 
semi-relativistic outflow. We simulate numerically the conversion of the MeV radiation into 
a vertical e ± outflow above a disc-like source. The outflowing e ± plasma becomes optically 
thick to Thomson scattering if the compactness of the 7-ray source exceeds ~ 30. The scat- 
tering by e then collimates the bulk of soft radiation along the disc axis, and the apparent 
bolometric luminosity of the disc depends strongly on its inclination to the line of sight. The 
anisotropic central emission may account for the lack of Fe Kcv lines in the X-ray spectra of 
bright radio-quiet quasars. The scattering in e outflows may also explain the orientation of 
optical polarization in non-blazar active galactic nuclei. 
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1 INTRODUCTION 

The X-ray activity of galactic black hole candidates (GBHs) and ac- 
tive galactic nuclei (AGNs) is usually attributed to accretion discs 
surrounding putative black holes. A special feature of accreting 
black holes is that their radiation spectra extend up to the 7-ray 
band. The 7-rays are thought to be produced in the innermost re- 
gion of the accretion disc (see, e.g., Svensson 1996 for a review). It 
may be a cold disc covered by an active corona in which particles 
are accelerated in magnetic reconnection events (Galeev, Rosner & 
Vaiana 1979). It may also be a hot two-temperature disc as pro- 
posed by Shapiro, Lightman & Eardley (1976). 

As soon as the radiation spectrum extends above the electron 
rest-mass energy, m e c? = 511 keV, the reaction 7 + 7 — > e~ +e + 
becomes inevitable and leads to production of an electron-positron 
atmosphere of the black hole (Herterich 1974). Guilbert, Fabian & 
Rees (1983) showed that if the 7-ray luminosity exceeds ~ 10~ 3 
of the Eddington limit for the black hole, the e atmosphere be- 
comes optically thick to Thomson scattering. The brightest accre- 
tion discs thus can be surrounded by optically thick e envelopes. 
The dynamics of an e envelope around a compact 7-ray source 
has been discussed in various contexts (e.g., Katz 1982; Zdziarski 
1984; Guilbert & Stepney 1985; Coppi & Lamb 1992; Illarionov 
& Krolik 1996; Li & Liang 1996), and it has been shown that the 
pairs should share the momentum of radiation and form an outflow 
from the source. 

In this paper, we study the formation of vertical e outflows 
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from disc-like 7-ray sources. Following the observed spectra of 
non-blazar AGNs and GBHs, we assume that the 7-rays contribute 
only a small fraction, Ss 0.1, to the bolometric luminosity, and the 
bulk of disc emission is in UV or X-ray bands. The created pairs 
then move in the dense soft radiation field which forces them to 
acquire an equilibrium velocity ~ c/2 such that the effective radia- 
tion pressure acting on the pairs vanishes (Gurevich & Rumyantsev 
1965). The equilibrium velocity is determined by the angular distri- 
bution of the radiation and can easily be calculated for an optically 
thin atmosphere. In the optically thick case, the problem becomes 
non-linear as scattering by moving pairs affects the angular distri- 
bution of the radiation above the disc. The corresponding transfer 
problem can be solved self-consistently to yield both the plasma 
velocity and radiation intensity as a function of optical depth in the 
e atmosphere (see Section 4). 

The density of the outflow is governed by pair injection due to 
7 — 7 interactions. We assume that the primary spectrum of 7-rays 
emerging from the source has a break at a few MeV due to strong 
absorption of more energetic photons inside the source. The bulk 
of 7 — 7 interactions outside the disc occur between MeV photons, 
just above the threshold. The 7-ray transfer is therefore non-linear 
as the 7 — 7 opacity at each height is determined by the MeV radi- 
ation field itself. To study the non-linear transformation of 7-rays 
into an e outflow, we have developed a code which calculates the 
hard X/7-ray transfer in a slab geometry. Besides the 7 — 7 reac- 
tion, the code approximately accounts for Compton scattering and 
annihilation emission of the created pairs, assuming that the pairs 
move with the typical velocity v = c/2. 

The main properties of the atmosphere, such as Thomson opti- 
cal depth and the power of e ± outflow, are determined basically by 
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one parameter of the source - the compactness parameter defined 
in Section 2.1. The results are weakly dependent on the specific in- 
ternal boundary conditions for the transfer problem. In our calcula- 
tions, the primary source is simulated in the simplest way: A homo- 
geneous disc emitting hard photons isotropically with a power-law 
spectrum around 511 keV. 

In Section 2, we discuss a general constraint imposed by 7 — 7 
absorption on the spectrum of an accretion disc. Then, in Section 

3, we calculate the rate of pair injection above the disk surface and 
discuss the formation of an optically thin e outflow. In Section 

4, we consider 7-ray sources with large compactness parameters, 
which are surrounded by optically thick e atmospheres, and cal- 
culate the e ± density and velocity profiles of the optically thick 
outflows. We also discuss pair acceleration at large distances from 
the source. The main conclusions are summarized in Section 5. 



2 GAMMA-GAMMA ABSORPTION 

7 — 7 absorption imposes a restriction on the maximum energy of 
the 7-rays emitted by a compact source (Guilbert et al. 1983). This 
provides a test for whether the observed 7-rays can be emitted by 
an accretion disc or if they should be associated with a relativistic 
jet or a more complicated geometry of the source (Illarionov & 
Krolik 1996). Such a test has been applied to 7-loud AGNs (e.g., 
McNaron-Brown et al. 1995; Becker & Kafatos 1995). 



2.1 The model 

Consider a homogeneous 7-ray emitting disc of radius R and as- 
sume that the hard X/7-ray emission emerging from the disc is iso- 
tropic and has a power-law spectrum around 511 keV with intensity 
(flux per unit energy per steradian) 



(1) 



where e = hv/m c c 2 is photon energy in units of the electron rest 
mass energy. The disc spectral luminosity is given by 

L s =L l£ - a , L x =2n 2 R 2 h. 

The main parameter of the model is the compactness parameter, 
h = (2) 



h = 2tt- 



m e c 3 i?' 

which can also be expressed in terms of the Eddington luminosity, 

Le = 2irr g m p c 3 /<tt, 
m p Li r g 
m c Le R ' 

where r g = 2GM/c 2 is the gravitational radius of a black hole of 
mass M. We assume the typical spectral index a > 1, then the 7- 
ray intensity steeply falls off towards high energies and h roughly 
corresponds to the total luminosity above 511 keV. This luminosity 
is a few per cent of the bolometric luminosity of the disc in UV and 
X-rays. 

2.2 Absorption inside and outside the source 

For a rough estimate of 7 — 7 absorption inside the source, take 
the 7 — 7 opacity of an isotropic radiation field (Gould & Schreder 
1967) 



^isotrop 



m c c z 



(3) 



where w(e) = 4irl £ /c is the spectral energy density of the radia- 
tion, and rj(a) is a numerical factor, e.g., r\ m 0.122, 0.072, 0.047 
for a=l, 1.5 , 2 respectively. The position of the spectral break due 
to absorption inside the source, e max , can be estimated from the 
condition K™° trop h ~ 1, where h is the disc thickness. Then one 
gets 



7T R 

2r]h ~h 



l/a 



(4) 



In the case h <JC R this is a much weaker constraint than that for 
a spherical source with h ~ R. However, equation (4) is not the 
final answer to the problem of 7 — 7 absorption because those 7- 
rays which have not been absorbed inside the disc-like source may 
be absorbed above the disc when passing through its radiation field 
(see also Zdziarski 1984). 

Let us evaluate absorption for a 7-photon of energy e emitted 
along the disc axis (z-axis) assuming for simplicity that the disc 
compactness is small (h < 1) and the non-linear effects are unim- 
portant. Then at height z the 7 — 7 opacity seen by our photon 
is caused by the radiation propagating freely from the disc within 
an angle 6> max = arctan(J?/z). Photons of energy e' streaming 
at angle 9 < # max can interact with our photon if e' exceeds the 
threshold 

2 

£thr = -f TT. (5) 

e(l — cosfl) 
The 7 — 7 opacity is then given by 



na T R J 



n / / — a— 1 
Q£ £ 



d(9sin<9 / 

x (1 — cos 9) a 11 {e, e 



(6) 



where <r 77 is the cross section for 7 — 7 pair production (Jauch & 
Rohrlich 1976) 



( \ — 3q " T 

J 77 (£c) — 
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(7) 



£ c — (e'/ £ thr) 1 ^ 2 is the energy of the interacting photons in the 
center-of-momentum frame. After integration we find 



= — 2 (l-coe< 



ct+2 a 



where 



77(a) 



a + 2 



y(£c) 



(8) 



<7T 



is the numerical factor already appearing in equation (3) and cal- 
culated analytically in Svensson (1987). Since the integral peaks 
at £ c ~ 1, the exact value of e m ax is not important and we set 
£max = 00. Now from equations (3) and (8) we get 

a + 2 



isotrop 



sin 



2a+4 



( flmax \ _ 1 / 1 

V^T) ~ 2F+ 2 \ 



Vz 2 + R 2 



This function is shown in Fig. 1 for a — 1 , 1 .5, 2. The opacity above 
the disc is strongly reduced because our photon moving along the 
disc axis does not encounter any radiation in the opposite direction 
which would dominate the 7 — 7 opacity in the isotropic case. 

It is interesting to compare the total optical depth seen by 
a 7-photon emitted along the disc axis, r 77 = J °° K 77 dz, with 
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Figure 1. The ratio of 7 — 7 opacity at height z above the disc to rop 
given by equation (3). The opacity is calculated for 7-rays propagating 
along the disc axis. R is the disc radius, a is the spectral index of the radi- 
ation. 

RK^° trop which can be (very roughly) associated with the 7 — 7 
optical depth of a spherical source of the same intensity, I E , and 
dimension, R. We get 

lT7 = sm 2 ° +4 (0/2) Ag 

7?4° trop J sm 2 9 

The dependence of this ratio on a is shown in Fig. 2. E.g., for 
a = 1, t 77 is « 30 times smaller than ii«*° trop . This result com- 
bined with equation (4) shows that the absorption constraint on 7- 
ray spectra is very sensitive to the source geometry and it can be 
much weaker than the usually used estimate for the spherical case. 

2.3 Dependence on inclination angle 

We now calculate the non-linear 7-ray transfer above the disc and 
study the dependence of the observed absorption break on the disc 
inclination to the line of sight. The bulk of 7 — 7 reactions oc- 
cur close to the disc surface (see Fig. 1), and we solve a simplified 
transfer problem in a one-dimensional approximation. In this ap- 
proximation, the radiation field is taken to be axisymmetric and 
uniform on any slice 2 = constant, and the effect of the finite 
size of the disc is treated for an optically thin atmosphere by trun- 
cating the radiation at angles 9 > 6* max (.z). We divide the slab 
z < R/2 above the disc into 25 horizontal layers, and assume that 
the radiation escapes freely from the last (upper) layer. The intrin- 
sic boundary condition is the isotropic power-law radiation given 
by equation (1). To get a stationary solution, I E (z,9), we take an 
initially empty atmosphere and let the 7-rays propagate from the 
disc. We then follow the evolution of the radiation field in the com- 
putational slab until a stationary state is established. Further details 
of the numerical method are given in Appendix A. 

Examples of spectra emerging at different angles from the disc 
with a modest compactness h = 25 and spectral index a = 1.5 are 



Figure 2. The ratio of the 7 — 7 optical depth seen by a gamma photon 
emitted along the disc axis to i?K 7 3 ° trop ; a is the spectral index of the 
radiation. 

shown in Fig. 3. We included only the 7 — 7 reaction in the transfer 
problem in order to see the efficiency of pure absorption. The other 
two non-linear effects, Compton scattering and annihilation emis- 
sion of the produced pairs, become important when h 30 as the 
atmosphere becomes optically thick, see Section 3.1. 

A qualitative dependence on the inclination angle 6 is seen in 
Fig. 3: Absorption is more efficient at large inclinations. At large 9, 
an emitted 7-photon travels a long path at small heights where the 
7 — 7 opacity is large. Such a photon encounters especially high 
opacity because the threshold condition (5) is weaker at large 9. As 
a result, the escape probability is strongly reduced. 



3 OPTICALLY THIN PAIR WIND 

Each 7 — 7 interaction produces one e pair. Therefore, the so- 
lution to the transfer problem in Section 2.3 automatically gives 
the pair injection rate, n^ 7 , as a function of the height z above the 
disc. The bulk of pairs are injected close to the disc surface and are 
immediately cooled by the soft radiation down to the Compton tem- 
perature of the radiation field, kTc ~ 1 — 10 keV. The time-scale 
for Compton cooling is given by 

tc — 3rn <=c _ 3 m e R / L \ -1 / R\ ^ 
8<jtw 16 m p c \Le/ \ r s J 

where w = L/ttR 2 c is the radiation energy density, L is the disc 
luminosity. Here only luminous discs are considered, and tc is 
shorter than the time-scale for pair escape, t csc ~ R/c. On the 
time-scale tc, pairs are also accelerated by the pressure of soft radi- 
ation. The e ± plasma is light and the gravitational force can be ne- 
glected as compared to the radiative force when L > (m e /m p )LE. 
Then pairs acquire an equilibrium bulk velocity for which the ra- 
diation pressure is balanced by the radiation drag. Near the disc 
surface, where the radiation field is semi-isotropic, the equilibrium 
velocity is about 0.45c (Icke 1989; see also Section 4.2). 
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£ = hi//m e c 2 



Figure 3. The 7 — 7 absorption observed at inclination angles 9 = 
0, 7r/6, it/3 from a disc with compactness £1 = 25 and spectral index 
a = 1.5. J| urface is the radiation intensity at the disc surface given by 
equation (1). I e is the intensity of the emerging radiation found by numeri- 
cal simulation of the 7-ray transfer above the disc. Compton scattering and 
annihilation emission in the e ± atmosphere are neglected. 



Table 1. The efficiency of 7-ray transformation into an e -wind 

h 5 10 15 20 25 30 

L±/Li 0.04 0.07 0.09 0.10 0.11 0.12 



3.1 Luminosity in pair rest mass 

In an optically thin wind, pairs escape before they can annihi- 
late since t csc is less than the time-scale for annihilation, t ann ~ 
(n+(7Tc) _1 . Then the flux of escaping pairs equals the column 
rate of pair production above the disc. The corresponding flux of 
e ± rest-mass energy equals 

r-R/2 



F± = 2m c c 



n 



11 



(z)dz. 



The total rest-mass energy carried away by the e wind can be 
estimated as 

L± = 2nR 2 F ± . 

In the limit of small compactness, 7-rays propagate almost 
freely in the atmosphere, and only a small fraction of Li is trans- 
formed into L±. This fraction is proportional to the 7 — 7 optical 
depth which, in turn, is proportional to l\ . Therefore, in the limit of 
small compactness L± /Li oc h ■ As h increases, 7 — 7 absorption 
significantly depletes MeV photons, which reduces the 7 — 7 opac- 
ity. The transfer becomes non-linear and L± jh\ saturates at Si 20 
%. Numerically calculated values of L±/L± for several compact- 
nesses h are given in Table 1. In the calculations, a spectral index 
a = 1.5 was assumed. 

Pairs escape with velocities ~ c/2, and the density of the at- 



mosphere can be estimated as n+ 
cal scattering optical depth of the e 
given by 

7T V Ll J 



TT 



- L ± /2-KR 2 m c c 3 . The typi- 
-cloud around the disc is then 



(10) 



Using Table 1, one can see that the atmosphere becomes optically 
thick when h ~ 30. 

In an optically thick wind, i ann would become comparable to 
t CBC . Then L± is limited by annihilation which depletes the particle 
population in the wind. The density of escaping pairs is controlled 
by the condition t ann ~ t CBC which gives n + ~ (er-r-R) -1 and 
the upper limit for L± (Phinney 1983; Guilbert & Stepney 1985; 
Ghisellini et al. 1992), 



Lmax 
± 



2-Km c c 3 R 



<7t 



(11) 



which corresponds to compactness l± ~ 27r. 



3.2 Annihilation emission 

Consider a plasma volume dV moving with velocity v = /3c. In 
the comoving frame, the number of annihilation photons emitted 
by this volume into solid angle dQ during time dt equals 

dN = — coTn+dVdidQ,, 
Ion 

where n+ = n+/7 is the pair density in the comoving frame, 
dV = 7dV, and 7 = (1 - p 2 y 1/2 . In lab frame, these pho- 
tons propagate within a solid angle dQ, — d£le~ 2 n , where e ann = 



account the invariance of four-volume, dVdt — dVdt, we get the 
annihilation power emitted in direction fl by a unit volume of a 
stationary wind, 



Q a nn(^) — 



e ann m e c dN 



3m e c 3 (TTn+ 



dVdtdO 167T7 5 (l-/3-ft) 3 ' 

Note that the fraction of photons emitted within a solid angle dfl 
equals df2£ ann /47r. The average energy of the annihilation photons 
equals ~frn c c 2 . 

A mildly relativistic bulk motion makes the wind emission 
considerably different from that of a cold source at rest: 

i) Q ann is significantly anisotropic even for a modest /3. E.g., 
for P = 0.5 the emission in the flow direction is 8 times larger than 
in the perpendicular direction. 

ii) The annihilation photons are blueshifted in the lab frame 
by a factor of [7(1 — /3 ■ O)] -1 . For a face-on-oriented observer 
and [3 — 0.5 the blueshift equals \/3. 

iii) The width of the annihilation line is determined by the 
velocity gradients in the wind rather than by the thermal motions. 
The bulk velocities are relativistic, so the line must be broad. 

The total power of the annihilation emission can be roughly 
estimated as 



\t~) * 

V k a nn / 



ttL±, 



where tt is the typical scattering optical depth of the atmosphere, 
see equation (10). In the limit of small compactnesses, L ann /ii 
falls off oc ll. Then annihilation of pairs produced inside the 
source is likely to dominate, and a narrow line at 5 1 1 keV might 
be expected. In the case of large compactnesses, where the atmo- 
sphere becomes optically thick, the line produced inside the source 
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must be scattered above the disc. The observed annihilation emis- 
sion is then dominated by the e outflow which produces a broad 
blueshifted line as discussed in next Section. 



4 OPTICALLY THICK OUTFLOW 

We now address the formation of an optically thick outflow by 
calculating the 7-ray transfer above a disc with a compactness 
ii > 30. The main conversion of 7-rays into e pairs occurs 
close to the disc surface, at heights 2 <C R (see Section 4.1). We 
therefore consider a simplified one-dimensional transfer problem, 
in which pair production is calculated in a slab < z < R/2. 
We thus neglect the boundary effects connected with the finite size 
of the disc. The radiation is assumed to escape freely at the outer 
boundary z = R/2. Further calculations verify that this is a rea- 
sonable first approximation. 

A magnetic field may affect the outflow dynamics if the mag- 
netic pressure exceeds the radiation pressure. The configuration of 
the accretion disc magnetosphere is, however, quite uncertain. We 
assume here that the magnetic field does not prevent the vertical e ± 
outflow. 

In an optically thick outflow t ann /t csc ~ {n+a^R)^ 1 < 1, 
i.e., pairs annihilate before they can escape. The bulk of e form 
a thermal population at the Compton temperature, kTc ~ 1 — 10 
keV, and only a small fraction form a non-thermal tail. The tail 
consists of those pairs which are braking after their creation to- 
wards thermal energies due to inverse Compton (IC) cooling by 
soft radiation and Coulomb collisions with the thermal e plasma. 
The stationary energy distribution of the braking particles peaks 
at semi-relativistic energies (Beloborodov & Illarionov 1995). The 
density of non- thermal pairs n+ th ~ nVtc is a small fraction, 
~ ic/iann, of the thermal population, n+. 

The pairs are created in the interactions between MeV pho- 
tons and they typically have initial Lorentz factors 7 ~ 2 — 3. 
IC photons produced in the slab therefore have modest energies 
£ic ~ 7 2 £s ^ 50 keV, where e s is the energy of a soft photon 
before upscattering, typically in the range 1 eV - 10 keV. (The non- 
thermal component is optically thin and repeated upscatterings are 
negligible.) It means that the kinetic energy of created pairs is just 
absorbed by the soft component of radiation and a negligibly small 
fraction is converted into hard radiation. We therefore neglect the 
presence of non-thermal pairs in the calculations of hard X/7-ray 
transfer and assume that the plasma is composed of the cold ther- 
mal pairs only. Then there are only two sources of hard photons: the 
primary source below the computational slab and the annihilation 
emission of pairs produced in the slab. 

In calculations of Compton scattering by thermal e we ne- 
glect their slow thermal motions and account for the fast bulk mo- 
tion which strongly affects the scattering process. The plasma has 
a bulk velocity near the equilibrium value v ~ 0.5c determined 
by the soft radiation which dominates the energy density of the at- 
mosphere. The equilibrium velocity can be found as a function of 
optical depth independently of the atmosphere density profile (see 
Section 4.2). The velocity varies from 0.3c to 0.7c. As a first step, 
we adopt = v/c = 0.5 in the simulations of the 7-ray transfer. 

4.1 Numerical model of outflow formation 

We have calculated several transfer models assuming an isotropic 
source (1) below the computational slab. Details of the code are 
given in Appendix A. 
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Figure 4. The density profile of the e ± outflow created above a disc with 
compactness l\ = 150 and spectral index a = 1.5 (solid curve). The 
dotted curve shows the density n* + which would correspond to the local 
balance between pair production and annihilation. 

We here consider the case where there is no flux from the 
7-ray source and set 71+ = at z = 0. Then the outflow is created 
above the source. An example of the density profile obtained for 
the atmosphere above a disc with a compactness h = 150 and a 
spectral index a — 1.5 is shown in Fig. 4. The equilibrium density, 
n*+{z), would correspond to local balance between pair production 
and annihilation, i.e., n+ = (3/8)(l — 2 )(JTCn+ 2 . The pair den- 
sity, n+ , strongly deviates from n*+ as a result of pair advection 
along the z-axis. The transfer equation in the stationary case 
can be written as 

dn+ 3 . 2x r * 2 21 

^"dT = 8^ ^ ^ T \- n+ ~ n+ J ' 

The first and second terms in the square braquets stand for e ± pro- 
duction and annihilation respectively. At large heights, advection 
results in that n+ > n+ and pair production becomes unimpor- 
tant: the outflow is dominated by pairs advected from the lower 
layers. The evolution of the pair density with height is then gov- 
erned mainly by annihilation. 

The steep decrease in n*+ shows that the rate of pair injection, 
h + J oc n+ 2 , steeply decreases with height as a result of strong 
degradation of the 7-raysQ. Roughly, the outflow is created in the 
"injection zone", z/R 0.2, where n+ < n+. The created pair 
plasma then outflows in the "annihilation zone", z <; 0.2, where 
n+ > 71+ . 

The intensity of the outgoing 7-rays is suppressed to such an 
extent that the equilibrium pair density becomes n+ ~ (<tt-R) -1 
which is also the typical density of the "false photosphere". At z ~ 

* Note that the radiation quickly becomes anisotropic in the low layers 
because of preferential depletion of 7-rays at large 8. The boundary effects 
are therefore reduced and the slab geometry may be a good approximation, 
in contrast to the optically thin case where the finiteness of the disc radius 
was important and we truncated the radiation at 8 > arctan(i?/z). 



© 0000 RAS, MNRAS 000, 000-000 



6 A. M. Belobowdov 



R we get n+ ~ (ctt-R) -1 , i.e., we are already in the transition 
zone between the optically thick e envelope and an optically thin 
wind from the envelope. At z > R the flow is no longer plane- 
parallel and the pair density is strongly reduced. Our code is not 
able to follow the transition to the optically thin wind. One can 
only estimate the flux of escaping pairs as ~ c/<7tR. The emerging 
luminosity in e rest mass is given by equation (11). 

We now briefly discuss the annihilation emission of the out- 
flow shown in Fig. 4. Only photons emitted perpendicular to the 
disc have a significant chance (~ 0.5) to escape the optically thick 
layers, and the bulk of photons emitted at large angles will be down- 
scattered or absorbed. The escaping photons will be detected by 
a distant observer as an annihilation line. In the numerical sim- 
ulation we obtained a line in the emerging spectrum located at 
£ann = [7(1 — ft cos 9)]~ where 3 = 0.5 is the atmospheric 
velocity and is inclination angle. With increasing inclination, the 
line flux as well as the continuum 7-ray emission decreases by 10 
- 30 times. The equivalent width of the line is ~ £ an n/2 at all in- 
clinations. A model taking into account the velocity gradient in the 
optically thick layers (see Section 4.2) would yield a broad annihi- 
lation bump of approximately the same equivalent width. 

4.2 Collimation of soft radiation 

The 7-ray flux converted into pairs above the disc is much smaller 
than the flux of soft radiation. The latter, being approximately equal 
to the total energy flux from the disc, is constant in the optically 
thick e outflow. However, scattering in the outflow changes the 
angular distribution of the soft radiation. 

The scattering mainly occurs in the optically thick layers, 
z < R, where the outflow is approximately plane-parallel. We 
consider the simplified one-dimensional problem of radiative trans- 
fer in a slab of vertically outflowing pair plasma atop an isotropic 
source of soft radiation. Each electron/positron moves in Compton 
equilibrium during its life-time in the e envelope. The equilibrium 
velocity at each height is determined by the local angular distribu- 
tion of the frequency-integrated intensity of the radiation, J. 

In the plane-parallel geometry, the intensity is a function of 
the angle 6 between the ray and the z-axis, and of the optical depth 
in the e ± envelope, dr = — 2n + aT<iz. The only parameter of the 
problem is the column density of the envelope, To /or . Here we ne- 
glect the polarization of radiation due to scattering (polarization is 
considered in Beloborodov 1998). The transfer equation then reads 

— — ^ - cos6»^ = (1 - Bcos6)(S - I), 

2n+a T cdt dr v ' A " 

< r < r , (12) 

The radiation scattered in direction fl is represented by the source 
function, 

^ ^ T (i-/3-») / J(n) ^ (n) 7 dn ' 03) 

where n is the direction of a photon before scattering and da /dfl 
is the differential cross section for Thomson scattering, e and e' are 
the photon energies before and after the scattering respectively, and 
they are related by the equation 



e'(l - /3 • O) = e(l - /3 ■ n), 



(14) 



which expresses the fact that the photon energy does not change in 
the comoving frame. 

The equilibrium velocity 8 is determined at each height by the 
equation 



(l + /3 2 )Ii-/3(/o+/ 2 )=0, 
where 



(15) 



Jm — 



47T 



7(61) cos 7 ™ 6 dSl. 



The equilibrium equation expresses the condition that the net flux 
of radiation in the comoving frame vanishes (see, e.g., Sikora & 
Wilson 1981), and in the lab frame it means that the radiation pres- 
sure is balanced by the radiation drag. Pairs keep their velocity near 
the equilibrium value which varies with height. The time-scale for 
relaxation to the equilibrium ~ tc is the shortest time-scale in the 
problem, therefore, the equilibrium velocity is a strong attractor in 
phase space, and we neglect deviations of /3 from that given by 
equation (15). The equations (12) and (15) then form a closed set 
of equations. 

We are looking for a stationary solution, 31/ dt = 0. The 
first and second moments of the radiation field are integrals of the 
problem, Ji(t) = constant and Izir) = constant, (in contrast 
to the classical problem of radiative transfer in an electron medium 
at rest, which has only h (r) = constant). The constancy of h 
can easily be checked by combining the first moment of equation 
(12) with the condition (15). J2(t) = constant is just another way 
to express the equilibrium condition: Anh/c equals the radiation 
pressure in the vertical direction, and its gradient vanishes in the 
equilibrium e ± envelope. 

Another special feature of the equilibrium radiative transfer is 
that the radiation density inside an e envelope cannot be strongly 
enhanced regardless the optical depth of the envelope. The trapped 
radiation is advected out of the optically thick layers with a veloc- 
ity ~ c/2, and the presence of the envelope around the disc weakly 
affects the radiation density, as can also be derived formally using 
h{r) = constant. The main impact of e ± is on the angular distri- 
bution of radiation. 

Here we adopt a simplified model for Thomson scattering. We 
will assume that the scattered radiation is isotropic in the comoving 
frame, i.e., a fraction AO./ An of the scattered photons have direc- 
tions within the solid angle dfl. The solid angle corresponding to 
dQ. in lab frame equals dQ. — 7 2 (1 — (3 ■ Cl) 2 d£l (see, e.g., Ry- 
bicki & Lightman 1979). The total scattering cross section equals 
a = <7t(1 — (3 ■ n), and the differential cross section for the scat- 
tering in the lab frame is given by 



<7t 



l-/3-n 



da_ _ 

dfl ~ 47T7 2 (l-/3-fl) 2 ' 
Equations (13) and (14) then give the source function 



S(0) 



Jo - 2(3h + B 2 h 



7 2 (1 -/3cos(9) 4 ' 

Using the equilibrium condition (15), we express h in terms of Jo, 
Ji, and /3, and rewrite S(6) as 

Jo - Bh 



S(8) 



(16) 



7 4 (1 -/3cos6») 4 

We solve the equations (12) and (15) with the source function 
(16) numerically. At the inner boundary (r = To), we set a source 
of isotropic soft radiation, which absorbs photons impinging back 
on the disc (the results are essentially the same if we assume re- 
flection at the inner boundary). At the outer boundary (r = 0), we 
assume free escape of the radiation. As initial conditions we as- 
sume that a semi-isotropic radiation is filling the atmosphere at all 
heights (curve 1 in Fig. 5b). This would be a stationary solution 
if the atmosphere were optically thin (the corresponding equilib- 
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Figure 5. The stationary solution for the soft radiation transfer in a plane-parallel e ± atmosphere with column density 3<r T : (a) Velocity profile. Parameter 
t is defined as dr = — 2n+cr T dz. The dotted line displays the equilibrium velocity in a semi-isotropic radiation field, (b) The angular dependence of the 
soft radiation intensity: 1 - as emitted by the disc, 2 - at the base of the atmosphere (r = 3), 3 - as emerging from the atmosphere (r = 0). The intensity is 
normalized so that = 1 for the radiation emitted by the disc. 



rium velocity c(4 — \/7)/3 w 0.45c is shown by the dotted line in 
Fig. 5a). The initial state evolves due to Thomson scattering in the 
atmosphere. The stationary solution J(r, 6), (3{t) that we are look- 
ing for depends only on the parameter To and does not depend on 
the density profile, n + (z). We therefore choose the simplest den- 
sity profile, n+{z) = constant. In the calculations, we use a grid 
which is homogeneous in the z- and ^-direction. The number of 
grid points is N z x N$ = 300 x 100. The chosen time step equals 
At = 0.3T O /2n + a T cN z . 

The resulting stationary solution for to = 3 is presented in 
Fig. 5a,b. In the lowest layers of the slab, the equilibrium velocity 
is less than 0.45c since, besides the semi-isotropic disc radiation, 
there is some flux in the opposite direction as a result of the ra- 
diation being backscattered by the outflow (see Fig. 5b, curve 2). 
In higher layers, the radiation gets collimated by scattering on the 
moving pairs, and the equilibrium velocity increases to ~ 0.7c. The 
angular distribution of the outgoing radiation is shown by curve 3 
in Fig. 5b. The distribution does not depend on To in the optically 
thick limit, and therefore does not depend on the compactness of 
the primary 7-ray source. The emerging radiation is significantly 
anisotropic, e.g., the intensity at 6 — tt/3 is ~ 7 times smaller than 
that at 9 = 0. 



4.3 Acceleration of escaping pairs 

The escaping pairs form an optically thin wind. At large heights, 
2 >• R, the angular size of the central luminous disc R ~ r g di- 
minishes and the escaping pairs are accelerated. Their equilibrium 
Lorentz factor is determined by the balance between pressure of 
the central radiation and drag by the radiation from the outer parts 
of the disc (Phinney 1987). With increasing height the density of 
the central radiation is reduced oc z~ 2 and redshifted in the co- 
moving frame as 7 _2 . The density of the diffuse radiation from 



the outer parts of the disc scales oc 2 ~ 3 , but it is blueshifted as 
7 2 . The resulting equilibrium Lorentz factor 7 cq ~ (z/r-g) 1 / 4 . At 
large distances from the black hole, the pair Lorentz factor falls 
below the equilibrium value because of radiation dilution, and sat- 
urates at 700 ~ (mpL/mcL^) 1 ^ , where L is the disc luminosity, 
Le is the Eddington limit. 

The collimation of the central radiation may increase 700 . The 
relevant collimation parameter, £, equals the ratio of the collimated 
intensity at 6 = to the intensity of isotropic radiation of the 
same net flux. In our case £ w 2.3, see Fig. 5b. Accounting for 
the collimation, the equilibrium Lorentz factor can be estimated as 
7cq ~ ( 1/4 (z/rg) 1//4 . At some height z*, 7 falls below 7 cq . At 
heights z » z„ the drag is negligible and the acceleration is deter- 
mined by the central radiation flux F = ^L/4tyz 2 , 



7 



(tt-F 



47 2 m c c 2 

Substituting 7 w c(d7/dz) and integrating this equation gives the 
asymptotic solution 

8 m e Le z 



z > z,. 



The falling of 7 below 7 oq and transition to the regime 7 <JC 7 eq oc- 
curs when the acceleration time-scale, 7/7, equals the escape time- 
scale, z/c. At z = z* one therefore has the condition 87 s (2/r g ) « 
(m p /me)((L/L-E). By matching the solutions for the two regimes 
7 = 7 eq and 7 <§C 7 oq at 2 = z, , one can estimate 2, and then 700, 



7o 



2/7 



1/7 



(17) 



The weak dependence on £ shows that 700 is not crucially affected 
by the collimation of the central radiation. To obtain large Lorentz 
factors ~ 10 — 20 observed in AGN jets, a collimation factor £ ~ 
10 2 would be required. 
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5 CONCLUSIONS 

MeV photons are strongly absorbed above a compact disk-like 
source and the observed 7-ray luminosity is suprressed, especially 
if the disc is viewed at large inclinations. The pair plasma created 
above the disc flows away in Compton equilibrium with the soft 
radiation field. The optical depth of the outflow depends on the 7- 
ray compactness, h, (see Section 2.1 for the definition of li). The 
outflow is optically thin if h < 30 and becomes optically thick if 
h i; 30. In the optically thin case, pairs escape without annihila- 
tion and the emerging luminosity in pair rest mass corresponds to 
the number of 7-photons absorbed above the disc. 

The optically thick outflow can be described as an envelope 
surrounding the disc where the created pairs have a life-time less 
than the escape time-scale. During their life-time pairs move in 
equilibrium with the soft radiation field with a velocity ~ 0.5c. 
Detailed calculations of the optically thick outflow show that its 
velocity increases from w 0.3c at the base to « 0.7c at the photo- 
sphere. In the photospheric layers, 2 ~ R, the pair density reduces 
to n+ ~ (ctt-R) -1 , R being the radius of the 7-ray emitting disc. 
A significant fraction of pairs created in these layers escape to pro- 
duce a luminosity in pair rest mass L± ~ 27rm c 3 i?/uT. 

The condition for formation of an optically thick e ± outflow 
is likely to be fulfilled in the brightest objects, such as quasars, 
especially if the black hole has a large spin. The bulk of the energy 
is then released in a very compact region, R w r g , and the outflow 
becomes optically thick if the luminosity above 511 keV exceeds 
a small fraction, ~ 3 x 10~ 3 , of the Eddington luminosity. The 
optically thick outflow will strongly affect the observed radiation: 

i) The luminosity of the central region becomes anisotropic 
as a result of scattering in the outflow. The apparent bolometric 
luminosity of the disc is then suppressed if the disc is viewed at 
large inclinations. 

ii) The outflow obscures the central region of the accretion 
disc and inhibits observation of reflection features, in particular the 
Fe Ka line. Reflected X-rays may be observed from the outer parts 
of the disc which is illuminated by the central source. However, be- 
cause of the strong anisotropy of the source, the reflection features 
should be weak. This may explain the lack of Fe Ka lines in bright 
radio-quiet quasars compared with lower luminosity Seyfert AGNs 
(Reeves et al. 1997). 

iii) The expected annihilation line is blueshifted and broad due 
to the velocity gradient in the outflow. The line increases the con- 
tinuum 7-ray flux by a factor of ~ 2 at 0.5 — 1 MeV, in contrast to 
the narrow annihilation feature at 5 1 1 keV usually expected from 
luminous 7-ray sources. 

Even an outflow of a modest optical depth, ~ 0.1, may affect 
the observed radiation by changing its polarization (Beloborodov 
1998): the polarization vector becomes parallel to the disc axis, 
in agreement with optical observations of non-blazar AGNs (An- 
tonucci 1992). 

Our calculations of the pair production around the accretion 
disc are limited to the one-dimensional slab approximation. Simu- 
lations of the outflow formation in cylindrical geometry may pro- 
vide new insights into compact 7-ray sources. The 7-ray emitting 
region of the disc has a finite radius and at large inclinations the 
slab approximation is not adequate as the bulk of observed hard ra- 
diation then comes from the boundary of the 7-ray emitting region. 
Simulations in cylindrical geometry would also allow the study of 
the interesting case where a source of GeV radiation is geometri- 
cally separated from the source of soft X-rays, being located in the 
centre of the surrounding X-ray ring (Illarionov & Krolik 1996). In 



such a geometry, the e outflow would be created in a cylindrical 
wall between the two sources. 
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APPENDIX A: NUMERICAL SIMULATION OF 
GAMMA-RAY TRANSFER 

The radiation field in a plane-parallel atmosphere is axisymmetric 
at any height z. The photon distribution over energies, e, and angles 
with respect to the z-axis, 9, is related to the radiation intensity by 



In the numerical simulations, we consider photon energies 0.1 < 
e < 10, and angles < 8 < n. The radiation field at each height is 
represented by photon numbers in N e x Ng = 50 x 20 cells which 
are logarithmically spaced in the e-direction and uniformly spaced 
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in the ^-direction. We checked that our results are not significantly 
affected by doubling the grid resolution. 

(i) 7 — 7 reaction 

To simulate the 7 — 7 process we need the corresponding effec- 
tive cross sections for the interaction between photons from cells 
(i s ,ie) and (i' e ,i' e ). Given an angle between the interacting pho- 
tons, £, and exact values of their energies, e and e', one can cal- 
culate the 7 — 7 cross section using equations (5) and (7) with 9 
replaced by £. To find the effective cross section for the interaction 
between two given cells we use the Monte-Carlo method. We cal- 
culate the cross section for randomly chosen values e, 9 and e', 6' 
within the cells (i e ,ie) and (i' e , i' e ) and randomly chosen azimuthal 
angles of the interacting photons. We then find the average cross 
section for 7 — 7 interaction between cells (i e ,i$) and (i' e ,i'g). Re- 
peating this procedure for each pair of cells, we get a matrix F of 
size N s x No x N e x No containing the effective cross sections. 
The rate of photon depletion in each cell due to 7 — 7 reaction is 
given by 

n 77 (i e ,i„) = cn(i e ,i e ) ^ F(i e , i e , i' e , i' e ) n{i' e ,i' e ). 



(ii) Compton scattering 

We also use the Monte-Carlo method to calculate the effective 
cross section for scattering from cell (i e ,ie) to cell (i' e , i' e ). In each 
Monte-Carlo event, we take a random photon energy e and angle 9 
within the cell (i e ,io), and calculate the total scattering cross sec- 
tion a = (1 — 3 cos 6)0. Here B is the plasma velocity (in units of 
the speed of light) directed along the z-axis, and a is the total cross 
section in the comoving frame (Jauch & Rohrlich 1976) 



3<7t 

~8f 



1 4 

2 + l 



2(1 + 2iy 



ln(l + 2e) 



where e = 7(1-/3 cos 9)e is photon energy in the comoving frame. 
The photon angle in the comoving frame is given by 

~ cos e - 3 



1-/3 cose' 



Then we perform a random scattering in the comoving frame. The 
differential cross section for Compton scattering is given by (Jauch 
& Rohrlich 1976) 



do" 3<tt 



dfi 



16tt [l + e(l-^)] 2 



1+M + 



l + e(l-At) 



where fi — ft • ft', ft and ft' being the photon directions before 
and after the scattering respectively. The photon energy after the 
scattering equals 



l+e(l-At)' 

Performing the inverse transformation to the lab frame, we get the 
energy, e', and angle, 9', of the scattered photon, and determine the 
corresponding cell (i' s , i' e ) into which it has been scattered. Repeat- 
ing 1.5 ■ 10 6 Monte-Carlo events for each cell (i £ .ig), we find the 
distribution of the scattered photons throughout all cells (i' e ,i'g) and 
get a matrix C(i e ,i$,i' e ,i' e ) containing the required effective cross 
sections. The photon sink and source due to Compton scattering, 
«c and fiQ respectively, are given by 

hc(i £ ,ig) = 2n+c n(i e ,ie)a(i e ,ie), 



nc(ie,ie) = 2n+c ^ C{i' e ,i' e ,i e ,ie) n(i' £ ,i' e ), 

where the bar over a denotes averaging within the cell (i e ,ie). The 
scattering can affect the total number of photons in the considered 
energy range 0.1 < e < 10 because scattered photons may end up 
outside this range. We checked that the code conserves the photon 
number, i.e., that the decrease of the total photon number due to 
scattering is equal to the photon flux through the boundary of the 
energy interval. 



(iii) Annihilation emission 

The probability that an annihilation photon has an angle within in- 
terval (0, 6 + dO) is given by (see Section 3.2) 



dP = 



sin 9d6 



2 7 2 (1 -Bcos8) 2 



Replacing d6 by ir/Ng and taking 9 = (n/Ne)(ie — 0.5), we 
have the probability P(io) for an annihilation photon to be in 
cell (i*,ie), where i* corresponds to the photon energy e ann = 
[7(l — /3 cos 6 1 )] -1 . The production rate of annihilation photons is 
then approximately given by 

n a nn(i e , ie) = 2n+ nn P(i e )c5 iei j , 

where n^ nn = (3/8)(l — B 2 )c<jtti\ is the rate of pair annihilation. 



(iv) The transfer equation 



The transfer equation reads 

dn 



dt 



n C + flann — C COS ( 



dz 



In the optically thin case (Sections 2.3 and 3.1), we neglect the 
terms and n ann , and we take into account the finite disc di- 
mension by truncating the radiation at angles 9 > arctan(i?/z). 
When simulating the optically thick atmosphere in Section 4.1, we 
account for all terms in the transfer equation and use the additional 
equation for e transfer, 



+ 



dn 
~~di 
where 



• 77 • ann n 

= n_j_' — n + — cp 



dn+ 



n 



= \y~ K J n- n (i E ,i e ). 



We model the atmosphere as a slab of height R/2 above the disc 
and divide it into 25 horizontal layers of equal thickness. Doubling 
the number of layers does not strongly affect the results. We set 
the boundary condition at z = as being a power-law isotropic 
radiation field. The outer boundary condition at z = R/2 is free 
escape of the radiation. The main parameter of the problem is the 
compactness, h, defined by equation (2). The other parameter is 
the spectral index of the primary radiation, a. To find a stationary 
state of the atmosphere we take an initially empty computational 
slab z < R/2 and follow the evolution of the radiation field and 
e ± plasma according to the transfer equation until a stationary state 
is established. The time step in the calculations equals 0.01i?/c. 
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